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We investigate the interaction of correlated electrons with acoustical phonons using the extended 
Hubbard-Holstein model in which both, the electron-phonon interaction and the on-site Coulomb 
repulsion are considered to be strong. The Lang-Firsov canonical transformation allows to obtain 
mobile polarons for which a new diagram technique and generalized Wick's theorem is used. This 
allows to handle the Coulomb repulsion between the electrons emerged into a sea of phonon fields 
{phonon clouds). The physics of emission and absorption of the collective phonon-field mode by the 
polarons is discussed in detail. Moreover, we have investigated the different behavior of optical and 
acoustical phonon clouds when propagating through the lattice. Initially the optical phonon fields are 
located at the lattice sites and do not spread out through the crystal and their evolution is limited by 
a time r, which means that the exchange of polarons by the phonon clouds is localized and occurs at 
one and the same lattice site. Hence, the renormalization of local polaron propagators due to optical 
phonons conserves the hopping matrix elements and the electronic band width. In the opposite case, 
when the phonon fields are of acoustical nature, the propagation of the fields is delocalized from 
the beginning and correlations lead band narrowing efi'ects. In addition, there is the possibility of 
electron transfer without being accompanied by phonon fields. In this case only the local electron 
propagators are changed while the tunneling matrix elements remain unaffected. In the strong- 
coupling limit of the electron-phonon interaction, and in the normal as well as in the superconducting 
phase, chronological thermodynamical averages of products of acoustical phonon-cloud operators can 
be expressed by one-cloud operator averages. While the normal one-cloud propagator has the form 
of a Lorentzian, the anomalous one is of Gaussian form and considerably smaller. Therefore, the 
anomalous electron Green's functions can be considered to be more important than corresponding 
polarons functions, i.e., pairing of electrons without phonon-clouds is easier to achieve than pairing 
of polarons with such clouds. 

PACS numbers: 78.30.Am, 74.72Dn, 75.30.Gw, 75.50.Ee 



I. INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tivity by Bednorz and Miiller^ the Hubbard model and 
related models such as RVB and t — J have widely been 
used to discuss the physical properties of the normal and 
superconducting states^^^. However, a unanimous ex- 
planation of the origin of the condensate in high tem- 
perature superconductors has not emerged so far. One 
of the unsolved questions is in how far can phonons be 
involved in the formation of the superconducting state. 
The aim of the present paper is to gain further insight 
into the mutual influence of strong on-site Coulomb re- 
pulsion using the single-band Hubbard-Holstein model^'* 
and a recently developed diagram technique^^-'^'^ . We 
consider now the most interesting case, namely super- 
conductivity of correlated electrons coupled to dispersive 
acoustical phonons. This investigation differs from our 
previous studies^''"^^ of electrons coupled to dispersion- 
less optical phonons, which was addressed by most other 
authors^""^^ . 

Because the interaction between electrons and phonons 
is strong, we include the Coulomb repulsion in the zero- 
order Hamiltonian and apply the canonical transforma- 
tion of Lang and Firsov^"* to eliminate the linear electron- 
phonon interaction. In the strong electron-phonon cou- 



pling limit, the resulting Hamiltonian of hopping po- 
larons (i.e., hopping electrons surrounded by phonon 
clouds) can lead to an attractive interaction among elec- 
trons meditated by the phonons. In this limit, the chem- 
ical potential, the on-site and inter-site Coulomb ener- 
gies as well as the frequency of the collective phonon- 
cloud mode (which is much larger than the bare acous- 
tical phonon frequencies) are strongly renormalized^^^^^. 
This affects the dynamical properties of the polarons and 
the character of the superconducting transition. We sug- 
gest that the resulting superconducting state with pola- 
ronic Cooper pairs is mediated by the exchange of phonon 
clouds and their collective mode during the hopping of 
the polarons. 



II. THEORETICAL APPROACH 

A. The Lang-Firsov transformation of the 
Hubbard-Holstein model 

The initial Hamiltonian of correlated electrons coupled 
to longitudinal acoustical phonons (the polarization in- 
dex is omitted) has the form 

H — He+ Hpf^ + He-ph, (1) 



where 



with 



He = ^ {t{j eo6ij} al^ttia- 

i ij 

^^P/^ = E'^k(&l6k+i), (3) 

k 

He-ph = ^gii- j)Qinj, (4) 

ij 
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Here ai„ (^L) ^-^^ annihilation (creation) operators of 
electrons at lattice site i with spin tr, 6k (^j^) are phonon 
operators with wave vector k; qi (pi) is the phonon co- 
ordinate (momentum) at site i, which is related to the 
phonon operators by 



Qi 



V2 



The Fourier representation of these quantities have the 
form 



bi 
Qi 
?k = 



-ikRi 
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(5) 



Pk = — (bl, - b-k 



In this Hamiltonian Uq and Vfj are the on-site and inter- 
site Coulomb interactions, t{i — j) is the nearest neigbor 
two-center transfer integral (which may be extended to 
include also next-nearest neighbor hopping of electrons) , 
g{i — j) is the matrix element of the electron-phonon in- 
teraction, So = eo — Ho, where sq is the local electron en- 
ergy and /io is the chemical potential of the system. Here 
and in the next part of the paper the Plank constant h 
is considered equal to one. The Fourier representation of 
t{j — i) is related to the tight-binding dispersion £(k) of 
the bare electrons with band width W, 



with R as nearest neighbor distance. Apparently the 
energy scale of the model Hamiltonian is fixed by the 
parameters W, U, g and Wk- The band filling n is an 
additional parameter. After applying the displacement 
transformation of Lang-Firsov^"', 



Hp = e^He- 



Ci(T — 6 di, 



^ia ^ "i. 



(6) 
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we obtain the polaron Hamiltonian in the form: 

Hp = Hp + Hpf^ + Hint, (8) 

H^ = Hi, Hl = eJ2 riia + Urn^m^, (9) 

Hint = ^tij - i)cl^Cia + ^Y'v^jmrij, (10) 



where 



A = at e"*'^' 



(11) 



= ^ ^ g(k)pke'''^* = J2Pj9i^j - R-')' (12) 



k k 

e = eo-ii, ii = iio + hVP'', U = Uo-VP\{13) 



and 



vpf^a-i) = 1 g(k)g(-k) ^-k(R,,-R,) 

^ iVV Wk 



(14) 



Hence, the effective intersite interaction is Va 



with Vi=j = 0. The frequency Wk of acoustical phonons 
is linear in k for sufficiently small wave vc;ctors. In or- 
der to have a reasonable expression for the parameter 
S'(k) of the canonical transformation, it is necessary that 
the condition g(k = 0) = is fulfilled. This condi- 
tion means that the movement of phonons with infinite 
wave length, which is equivalent to the macroscopic dis- 
placement of the system, cannot influence its properties 
and must be omitted. Therefore, the Fourier representa- 
tion of the direct attraction mediated by phonons must 
also vanish in this limit: yp'*(k = 0) = 0. It IS im- 
portant to note that the Fourier representation of the 
Coulomb part of the inter-site interaction must also van- 
ish for vanishing wave vector: F'^(k = 0) = as a 
consequence of required charge neutrality of the system. 
Therefore, the resulting direct interaction between elec- 
trons, ]/(R) = y'=(R) - yP'*(R), fulfills y(k = 0) = 0. 
This will be used when analyzing the corresponding dia- 
grammatic contribution. 

When deriving the polaron Hamiltonian, it was nec- 
essary to include the shift of the polaron coordinate gk 
by 



e^gke ^ = gk - 



1 



^ k 



g{\s)nie 



ikRj 



which helps to eliminate the linear electron-phonon in- 
teraction. 
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The polaron Hamiltonian is by nature a polaron- 
phonon operator because the new creation and annihi- 
lation operators Cj^ and Cio- entering Hp must be inter- 
preted as operators of polarons, i.e., electrons dressed 
with displacements of ions that couple dynamically to 
the momentum of acoustical phonons. In the zero-order 
approximation (omitting Hint) polarons are localized and 
phonons are free with a strongly renormalized chemical 
potential // and on-site interaction U. This last quantity 
can become negative if the phonon mediated attraction 
V^^ is strong enough to overcome the direct Coulomb re- 
pulsion. The first term of the perturbation operator Hj^i 
describes tunneling of polarons between lattice sites, i.e., 
tunneling of electrons surrounded by clouds of phonons. 
The second term of this operator describes the renormal- 
ized polaron-polaron inter-site interactions. 



B. Averages of electron and phonon operators 

One problem is to deal properly with the impact of 
electronic correlations on the polaron formation involving 
operators like (11) for the electron and phonon subsys- 
tems. This can be done best by using Green's fimctions 
provided one finds a way to deal with the spin and charge 
degrees of freedom. In order to achieve this in the limit 
of large U, the Hubbard term can be included in the 
zero-order Hamiltonian. As a consequence, conventional 
perturbation theory of quantum statical mechanics is not 
an adequate tool because it relies on the expansion of the 
partition function around the noninteracting state using 
Wick's theorem and conventional Fcynman diagrams. 

In order to have a systematic description of corre- 
lated electrons, Hubbard^ proposed a graphical expan- 
sion around the atomic limit in powers of hopping in- 
tegrals. This diagrammatic approach was reformulated 
by Slobodyan and Stasyuk^^ for the single-band Hub- 
bard model and independently by Zaitsev^^ and further 
developed by Izyumov^'^. In these approaches, the com- 
plicated algebraic structure of the projection or Hubbard 
operators was used. 

We have found an alternative way in the sense that our 
diagram technique involves simpler creation and annihi- 
lation operators for electrons at all intermediate stages 
of the theory and Hubbard operators only when evalu- 
ating final expressions^"^''. In this approach, averages of 
chronological products of interactions are reduced to n- 
particlc Matsubara Green's functions of the atomic sys- 
tem. These functions can be factorized into independent 
local averages using a generalization of Wick's theorem 
(GWT) which takes strong local correlators into account, 
see Refs. 9,10 and 17 for details. Application of the GWT 
yields new irreducible on-site many-particle Green's func- 
tions or Kubo cumulants. These new functions contain 
all local spin and charge fluctuations. A similar linked- 
cluster expansion for the Hubbard model around the 
atomic limit was recently formiilated by Metzner^^. But 
in the latter work the Dyson equation for the renormal- 



ized one-particle Green's function was not derived, nor 
the correlation function which appears as main element 
of this equation. It is the purpose of this paper to check in 
how far we can use the GWT for the extended Hubbard- 
Holstein model given by Eq. (1). 

With respect to the transformed Hubbard-Holstein 
model, phonon operators are averaged using ordinary 
Wick's theorem by taking into account the factorization 
of the phonon partition function in k space of phonon 
wave vectors. We define the temperature Green's func- 
tion for the polarons in the interaction representation by 

Gp(x,a,T|x',a',T') = - (T Cx.(t)cxV' (^')C^(/3))o , (15) 
with 

cx.(r) = e""^c^^e-""\ Cx.(r) = e^^^c^.e-^''^ 

for the polaron and phonon operators, respectively, with 
= Hp + H^f^. Instead of i, j we now use x, x' as 
site indices; r, r' are imaginary time variables with < 
r < /?; T is the time ordering operator and /? the inverse 
temperature. The evolution operator is given by 



U{l3) = T exp 



■/3 



drH 



int\ 



(16) 



The statistical averages (...)o are evaluated with respect 
to the zero-order density matrix of the grand canon- 
ical ensemble of localized polarons and free acoustical 
phonons. 



n 



n 



-/3uikbl_bk 



(17) 



The subscript c in Eq. (15) indicates that only connected 
diagrams have to be taken into account. The polaron 
part of the density matrix (17) is factorized with respect 
to the lattice sites. The on-site polaron Hamiltonian con- 
tains the polaron-polaron interaction which is propor- 
tional to the renormalized parameter U. Therefore, this 
Hamiltonian can be diagonalized only by using Hubbard 
operators^. At this stage no special assumption is made 
about the value of the quantity U and its sign. So we 
can set up the eqiiations of motion for the dynamical 
quantities in this general case. Wick's theorem of weakly 
coupled quantum field theory can be used to evaluate 
statistical averages of phonon operators, including, the 
propagator of phonon clouds. 



C. Phonon-cloud propagators 

The zero-order one-phonon Matsubara Green's func- 
tion has the form 

a{x,x') = a(x-x'lT-T') = (T7rx(r)7rx'(T'))o 



4 



= ^El5(k)l'cosk(x-x') 

k 



cosha;k(/3/2 



r'l) 



sinha;k/3/2 



(18) 



with 



7rx(T) 



j 



Pj(T)5(R-i -Rx). 



Here X is again the position and r the imaginary time 
while X in Eq. (18) stands for (x, r) 

This function makes an essential contribution for small 
values of distances |x — x'| and |t — t'| close to zero or /?. 
For X = x' the minimum value of this function is obtained 
for |t — r'l = /3/2. Since all approximations in this paper 
concern the strong-coupling limit of the electron-phonon 
interaction, we will use the series expansion of a{x,x') 
near r = and t = (3: 



a(0|T) 



(t(0|0) - uJcT, r > 

a(0|0)+We(r-/?), r<(3 



with 



2N 



(19) 



(20) 



as collective phonon cloud frcquency^"^'^^. Besides the 
one-phonon propagator we have also many-phonon cloud 
propagators. There are two kind of one-cloud propaga- 
tors, of which (f>{x\x') is the normal-state one and ip{x\x') 
the anomalous one of the superconducting state, given by 



^(x|x') 


= -/-(x 


-x'lr-r') = 






= exp 


(-i(T[^x(r)- 


-^x'(r')]')o) 




= exp[ 


-a{0\0) + a{x 


-x'|r-T')], (21) 


ip{x\x') 


= <y5(x 


— x'|r — r') = 


^'Pgii-=c(-r)-|-i7r^/(r')^^ 




= exp 


(-^(Tkx(T)^ 


-'rx'(r')]')o) 




= exp[ 


-cr(0|0) - C7(x 


-x'|r-r')]. (22) 



For the first function the maximum value of the one- 
phonon propagator (t(x|t) is favored while for the second 
one the corresponding minimum value is preferred. The 
Fourier representations in r-space have the form 



P n 

*""-^(iO„), 



n 



where 



J o 
J o 



^^giQrg-<7(0|0)-<7(0|r)_ 



(23a) 
(23b) 

(24a) 
(24b) 



Here is r2„ the even Matsubara frequency n„ = 27rn//3. 
In order to find the Fourier representations of these func- 
tions we have used the peculiarities of the cr-propagator 
in the strong-coupling limit of the electron-phonon inter- 
action. As proven in Appendix A, the first propagator 
can be written as 



0(x|t') ~ (/.(x)<^(r), </)(x) « 5: 



x,0 



with 



(q) 



(25) 



A more realistic value for 0(q) is obtained by using the 
dependence of (7(x|r) on small values of x. In this more 
precise approximation we find 



<^(q) 



where 



27r\ 
^1/ 



3/2 



-qV(2<7i) 



<^(x) 



-aixV2 



(26) 



CTi 



^El^(k)l'k^cothia;k/?. (27) 



This result has been obtained by an expansion of cos kx 

in terms of x. We also assume that ^(k) depends on 
k only through its modulo |k|. Then the Fourier rep- 
resentation of the normal phonon cloud propagator is 
a Lorentzian and therefore the time dependence of this 
phonon cloud corresponds to that of an oscillator with 
the large collective frequency Uc- For the anomalous one- 
cloud propagator (p{x\x') we obtain in this approximation 
a Gaussian representation, see Appendix A: 



= -\/27r/cr2 exp [^pCln 

- am - <7(0|/J/2) - (0„)V(2a2)] , (28) 



where 



a2 = a" (01/3/2). 



(29) 



The space dependence of (p(x|zri„) is more complicated 
compared to the space dependence of (f>(x\iCl„) because 
we cannot restrict the discussion to small values of |x|. 
In the following we will discuss many-cloud propagators, 
both in the normal and superconducting states. We start 
with the two-cloud propagators [as before, x = (x, r)]: 
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(f>2{xi,X2\x3,X4) = (T exp[i [7rxi(Ti)] +7rx2(r2) - 7rx3(T3) - 7rx4(r4)])o 

= exp (-^(T [tTxiCti) + 7rx2(T2) - 7rx3(r3) - 7rx4(T4)]^)o) 

= exp (S(xi,Ti;x2,r2|x3,T3;x4,r4)), (30) 

iP2{xi,X2,Xs\x4) = (T exp[i[7rxi(Ti)] +7rx2(T2) +7rx3(T3) -7rx4(r4)])o 

= exp (-5(T [7rxi(ri) + tt^^{t2) + 7rx3(r3) - 7rx4(T4)]^)o) 

= exp (S(xi,Ti;x2,r2;x3,T3|x4,r4)), (31) 



where 



S(xi,Ti;X2,r2|x3,T3;X4,r4) = (t{xi - XsUn - T3I) +<t(xi - X4||ti - T4\) 

+ £7(x2 - X4||r2 - T4I) +a-(x2 - X3IIT2 - rsl) - £7(xi - X2||ri - T2I) - a{x3 - X4\\t3 - T4\) - 2£7(0|0), (32) 

S(xi,ri;X2,T2;X3,r3|x4,T4) = ct(xi - X4||ti - T4I) + Cr(x2 - X4||t2 - T4I) 

+ (t(x3 - X4||-r^ - T4I) - cr(xi - X2IIT1 - T2I) - cr(xi - xslln - T3I) - ct(x2 - X3IIT2 - T3I) - 2(7(0|0). (33) 



The following relations exist between two- and one-cloud Green's functions: 








<?i>2(a;i,a;2|a;3,a;4) = (/)(xi X3)(/)(x2 X4) exp 


[(T(a;i X4) + cr(x2 X3) 


- a{xi\x2) 


- cr(x3 X4)] 


(34a) 


= (I>{xi\x4)(p{x2\x3) exp 


[cr(xi|x3) + a{x2\x4) 


- a{xi\x2) 


- a{X3\X4)] , 


(34b) 


</?2(a:;i,a;2,a;3 2:4) = ^(xi a;2)0(a;3 ^4) exp 


[a{xi\x4) + cr{x2\x4) 


- (7(X1 X3) 


- 0-(x2 X3)] 


(35a) 


= (^(aJi a;3)^(a;2 a;4)exp 


[a{xi\x4) +a{x3\x4) 


- Ct{xi\X2) 


- a-(X2|X3)] 


(35b) 


= (^(a;2|x3)0(a;i|a;4)exp 


[(t{X2\X4) +a{X3\X4) 


- a{xi\x2) 


- <t(xi|X3)] . 


(35c) 



Many-cloud phonon propagators will be present in all di- 
agrams of the thermodynamical perturbation theory to 
be formulated here. As above equations show, all sites 
of the diagrams are joint and appear to be connected 
in the presence of acoustical phonons. In order to clas- 
sify the diagrams as connected and disconnected ones, it 
is necessary to have the analogy of Wick's theorem for 
many-cloud propagators similar to the theorem we had 
formulated for correlated electrons^'^"'^''. In the absence 
of such a theorem we cannot prove the existence of a 
linkcd-clustcr theorem for the thermodynamical poten- 
tial and for other extensive quantities. 

This problem has been discussed in detail in Ref.^°, 
however, only now we arc able to present a solution. In 
order to obtain this solution, we observe that the two- 
cloud functions determined by Eqs. (34) and (35) have 
their maximum values when the arguments of the normal 
one-cloud functions 4>{x\x') coincide (x = x') and the cor- 
responding exponential factors close to these arguments 
approach one. There are several possibilities to achieve 
this and all of them have to be taken into account. We 
assume that as main approximation the following expres- 
sions for the two-cloud propagators will result, 

02(X1,X2|X3,X4) = (/'(Xi|x3)(/)(X2|X4) 



-I- (^(Xi|x4)(/)(X2|X3) 

-I- 02'^(X1,X2|X3,X4), (36) 

V2(X1,X2,X3|X4) = V'(a^l|X2)^(X3|X4) 

-I- (/9(Xi|x3)(;!)(x2|X4) 

+ (/?(X2|X3)0(X1|X4) 

-I- Qi2''(xi,X2,X3|x4). (37) 

These last equations also define the irreducible parts of 
the two-cloud propagators or phonon-cloud cumulants. 
In the strong-coupling limit the irreducible functions are 
small and can be omitted as shown below. The validity 
of this statement is discussed in Appendix A, in which 
the Fourier representation of the normal two-cloud prop- 
agator, 

02(xi,zili;x2,if22|x3,if23;x4,if24) = / ... / rfri...rfr4 

Jo Jo 

X exp (ifiiTi -|- ifl2T2 — ii^3T3 — ifl4T4) 

X(/)2(xi,Ti;x2,r2|x3,T3;x4,r4), (38) 
has been calculated in the strong-coupling limit leading 
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to 

(j!)2 (xi , ifii ; X2 , iSl2 1 X3 , i ^3 ; 2:4 , ^^^4) 

~ 0(xi - ■X.2\iill) fciO;, '/'(X2 - X4|i02) (50204 

+ 0(xi - X4|zfii) (5oi04 <i'(x2 - X3|jrj2) (5o203-(39) 

The last equation shows that in this hmit the irreducible 
function is not relevant and Wick's theorem has a simple 
form, which does no contain significant irreducible con- 
tributions. Similarly we obtain for ip2 a form without 
irreducible contributions, 

<^2(xi, iOi; X2, iCt2; X3, i^lalxi, i^li) 

— I I ^j-^ ^—{iQlTl+iQ2T2+iQ3T3—i^4,Ti) 

Jo Jo 

X (/32 (xi , n ; X2 , T2 ; X3 , rs |x4 , T4 ) (40) 

~ ip{xi - X2|iril) (5o2,-Jli0(x3 - X4|iri3) (^03,04 

+ (/?(xi - X3|iOi) (5o3,-ni <^(x2 - X4|iri2) (5o2,04 

+ (f{x2 - X3|if^2) foa -Ha " y^4\i^4) fci,04- (41) 

These results correspond to our preliminary estimates 
that the irreducible parts in Eqs. (36) and (37) can be 
omitted because they arc not important in the strong- 
coupling limit, see Appendix A. Hence, without the ir- 
reducible parts the equations assume a form correspond- 
ing to Wick's theorem applied to two-cloud propagators. 
This can easily be generalized to the case of a larger 
number of clouds. Thus, there is an analogy of having 
a generalized Wicks's theorem for the case of correlated 
electrons^' and a corresponding theorem for correlated 
phonon clouds. This allows us now to develop a thermo- 
dynamical perturbation theory for correlated electrons 
interacting strongly with phonons. 

As is shown below the tunneling of polarons between 
lattice sites can be accompanied by either preserving or 
by exchanging phonon clouds. In the strong-coupling 
limit these clouds are heavy, therefore, in case of pre- 
serving the cloud, the effective matrix transfer matrix 
element is considerably diminished leading to band nar- 
rowing effects. In the other case, when clouds are ex- 
changed, the transfer matrix element and the electronic 
band width remain unchanged. 



III. POLARON GREEN'S FUNCTIONS 



A. Local approximation 

The zero-order one-polaron Green's function is given 

by 

GO(ar,ar') = -(Tcxa(r)cxV'(r))o 

= -(Taxa(T)axV'(T))o H^, r|x', t') 

= ^x,x'<5a,.'G°(T-r')<^(r-r'), (42) 

where x stands now for x = (x. a. r). In order to discuss 
the influence of the collective mode on Gp{x, x'), we write 



down its Fourier transformation by making use of Eqs. 
(25) (see Ref. 19): 



1 

Z~o 



dTe'^"^G%{T), (43) 



Ef) — E„ — LUc 



+ 



+ 



NjiOc) (e"^^"+e-^^-) 
e-/3-E-. +N{tj^) (e-Z^-E- + e-'5-E2) 



E. 



E2 



-0E2 



+ 



-I3E2 



E. 



E2+0JC 

where w„ is the odd Matsubara frequency and 

Zo = 1 + e-f"^- + e-^^-- + e-^''^ 
Eo = 0, E±^ = e, E2 = U + 2e, 
n{e) = {ef'' + N{lu^) = (e^'^'^ ^ 



(44) 



(45a) 
(45b) 
(45c) 



Equation (44) shows that the on-site transition energies 
of polarons are changed by the energy of the collective 
mode. The delocalization of polarons due to hopping and 
intersite Coulomb interaction leads to broadening of the 
polaronic energy levels. The polaron propagator has the 
following antisymmetry property: 

G%{iu)n; e; We) = -G° (-iw„; -s - U; Wc) (46) 



B. Expansion around the atomic limit 

We will now investigate polaron delocalization under 
the influence of Hint, in Eq. (10) by making use of thcr- 
modynamical perturbation theory in the interaction rep- 
resentation. The averages of chronological products of in- 
teractions are reduced to n-particle Green's functions of 
the atomic system, which can be factorized into indepen- 
dent local averages of electron operators and chronolog- 
ical products of phonon operators. The procedure relies 
on a generalized Wick's theorem for electron operators, 
which takes into account the strong local electronic corre- 
lations, and Wick's theorem for phonon cloud operators. 
In addition to the normal one-polaron propagator in Eq. 
(15), we will also investigate the anomalous propagators 
defined by 

Fpix\x') = -(Tc,c,.[/(/3))S, (47a) 
Fpix\x') = -(Tc,c,.t/(/3).)S (47b) 

As before, x stands for (x, a, r). The easiest way to estab- 
lish (47) is to make use of a local source term of Cooper 
pairs. 
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which is added to the local Hamiltonian (2) and switched 
off at the end of the calculation. 

In first order perturbation theory the contributions to 

the normal polaron Green's function (15) and anomalous 
Green's function (47a) are shown in Fig. 1(a) and Fig. 
1(b), respectively. 

The diagrammatic elements are self-explanatory (see 
caption of Fig. 1). Since the correlation functions 
r2(a;, x'\i) and ^^(x, x'\i) and the two-particle irreducible 
fimction *'^[a;i, X2|a;3,a;4] are local quantities, all site 
indices are equal: 

T°{x,x'\i) = ^x,x'^x,irO(a,r;a',r'|ri), (48a) 
'^Ux,x'\i) = 5x,x'5x,i9^(a,r;a',r'|ri), (48b) 

where 

rO(a,r;a',r'|ri) = {T a^{T)a„,{T')n{n))o 

- (Ta^(T)a.'(r'))o(n(ri))o, (49a) 

Q2{o',T;a' ,t'\ti) = {T acr{T)a„^{T')n{Ti))a 

-{Ta^{T)a„,{T'))o{n{Ti))o. (49b) 

These functions can be compared with the two-particle 
irreducible quantities of Ref. 9,10 defined by 

G2^''[X1,X2\X3X4\ = (Jxi,X2^xi,x3^xi,X4 

X G^^'^[aiT-i;(T2T2\(T3T3;(74T4\, (50a) 
G^''''[xi,X2,X3\X4\ = ^xi,X2'5xi,X3^xi,X4 

X G^'''[<Tiri;£72T2;<73T3|£74r4], (50b) 

where 

<J2"' [criTi;a2T2\a3T3; aiTi] 

= {Ta^^ {Ti)a^2 {T2)aa3 (^3)0^4 i'^4:))o 

- (Ta^i (Ti)acr4(T4))o {T^ a„^{T2)aa3{T3))o 

+ (Taai(ri)a(T3(T3))o (Ta<^2(T2)aa4(T4))o- (51) 

The second function, G2^^[xi,X2,X3\x4\, can be ob- 
tained from Eq. (51) by replacing the operator 0,73 (T3) 
by aa.^{T3) and corresponds to an anomalous supercon- 
duc'ting contribution. Between the non-full cumulant 
T%{a,T\a' ,t'\tx) and the full one exists the follow- 
ing relation 

rO(a,r;a',r'|ri) = - GO-[aT; c7iTi|<7ir+; aV] 

+ ^G"(aT|airi)G"(airi|a'T').(52) 

In second order perturbation theory we have to deal with 



more complicated functions like 

T3{x,x'\ii,i2) = ^x,x' S^,h ^x,i2 r3(crr;a-'T'|Ti,T2)), 

(53a) 

r3(ar;aV|Ti,T2) = (T a<,(T)a^/ (T')n(Ti)n(T2))o 

- (Taa(r)a<^/(T')n(ri))o («(t2))o 

- (Taa(r)a<,/(T')ra(r2))o (n(ri))o 

- {aa{r)aa'iT'))Q (n(Ti)n(T2))o 

+ 2{aa{T)aa' {T'))o{n{Ti))o{n{T2))o, 

(53b) 

with the following relation between the non-full and full 
three-particle cumulants, 

T3{aT;a'T'\Ti,T2) 

= - 5Z G3''''[aTaiTi;a2T2\aiT^;a2T} icr'r'] 

(T1(T2 

- ^ G\aT\a,n)G\a,n\a2T2)G\a2T2\a'T') 

- J2 G°(fTT|a2T2)G0(a2T2|aiTl)G0(<7lTl|<7V'),(54) 

where — ti + and T2 = T2 -I- 0. The diagrams of 
Fig. 1 contain also this element which is marked as a cir- 
cle at the end of the long-dashed lines. This element is 
the average of the electron number operator n. The con- 
tribution of such diagrams is proportional to y(k = 0) 
which is equal to zero. Hence, they can be omitted. In 
still higher order perturbation theory these circles appear 
with TO long-dashed lines ending on them, which repre- 
sents the TO-order Kubo cumulant for the mth degree of 
the number operator n: 

n^^ = {{h-{h)f)-mn-{n)f)f. (55) 

The first diagram of the right-hand part of Fig. 1(a) is 
a local normal electron propagator G^{x\x') represented 
by a solid line, which is renormalized by the phonon- 
cloud propagator 4){x\x'). The corresponding first dia- 
gram in Fig. 1 (b) is the anomalous electron Green's func- 
tion renormalized by the phonon-cloud. In this case the 
renormalization is determined by the anomalous function 
ip{x\x'), which is smaller than the corresponding contri- 
bution from the propagator (j){x\x'). The diagrams (b)- 
(e) of Fig. 1(a) and (b)-(g) of Fig. 1(b) take into account 
intersite tunneling, of which the corresponding last dia- 
grams describe direct polaron-polaron intersite interac- 
tion. 

Among these diagrams we have also weakly connected 
ones, e.g., diagram (b) in Fig. 1(a) and diagrams (b) and 
(e) in Fig. 1(b), which can be divided into two parts 
by cutting one tunneling line. All other diagrams are 
strongly connected. All contributions of Fig. 1(b) for the 
anomalous Green's functions contain the phonon-cloud 
propagator Lp{x\x') and are therefore less important than 
the corresponding diagrams of Fig. 1(a). 
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(h) 

FIG. 1: The simplest diagrams contributing to the normal (a) and anomalous (b) one-polaron Green's functions. Solid lines with 
arrows in same direction represent normal (G") and lines with arrows in opposite directions anomalous {F^\ F") propagators, 
respectively. Short-dashed lines are for the hopping matrix elements t{i — j), long-dashed lines represent the direct polaron- 
polaron interactions V{i — j), the wiggly lines stand for the normal phonon (cloud) propagators (l){x\x'), and the zigzag lines 
are the anomalous phonon (cloud) propagators ip(x\x'). Rectangles depict the correlation functions r2 and 



In second order perturbation theory there are both 
groups of weakly and strongly-connected diagrams. The 
weakly connected ones are the simple repetition of the 
first-order diagrams. Therefore, in Figs. 2-4 only the 
strongly-connected diagrams for the normal and anoma- 
lous polaron Green's functions are shown. 

The diagrams (a)-(i) in Figs. 2 and 3 originate from 
polaron tunneling processes. In the group of more com- 



plex diagrams the contributions from (f) appear only be- 
cause of the electron-phonon interaction. Diagrams (g)- 
(k), each with a rectangle, are special contributions char- 
acteristic of strongly correlated electron systems. The 
diagrams (h) and (i) contain the phonon (cloud) cor- 
relation function (f>3{xiii2\jij2x'), each consisting of six 
terms with different order of phonon propagation. To 
illustrate this, the corresponding diagrammatical repre- 



G° ji ii G° 72 h G" 

(b) 




(e) 



FIG. 2: Strongly-connected diagrams of second order perturbation for the normal polaron Green's function. 



sentation for one of the anomalous Green's functions is 
shown in Fig. 4 (without resolving the structure of the 
less important polaron-cloud propagator ipn)- 

With respect to the full normal Green's function Gp we 
now introduce the correlation function Zp(x\x') in anal- 
ogy with previous works^"^^. In diagrammatical form 
typical contributions to this function are given by (c)-(f) 
in Fig. 1(a) and by all contributions from Figs. 2 and 3 
and by corresponding strongly-connected contributions 
from higher order perturbation theory. If we then con- 
nect the strongly-connected diagrams in all possible ways 
by the weakly connected diagrams with non-renormalized 
tunneling matrix elements, we are able to formulate a 
Dyson-like equation, 

Gpo.o.'(x,t|x',t') = Apcra'(x, r|x',r') 

+ / ^"^1 X! X! '^1^1' "^i) 

XiX2 (Ti 

X t(xi - X2) Gpfficr'(x2,Tl|x',T'), (56) 

where 



Ap(a;|a;') 



G^(x|x') 



which in Fourier representation, 



has a simple algebraical form: 



Gpaa{k\iLUn) 



Ap^^{k\iLUn) 



1 ~ e(k) Ap^a{k\iujn) ' 



(59) 



Zp{x\x'), 



(57) 



with the tight-binding dispersion e(k) of the bare elec- 
trons defined before. 

Here we have assumed a paramagnetic ground state 
and spin conservation a' — a. It is important to note 
that the form of the Dyson equation is the same for both 
the superconducting and normal states of the system. 
The states differ with respect to the correlation func- 
tion Zp. This situation is somewhat different from the 
Hubbard-Holstein model for optical phonons^^, where in 
the superconducting state normal and anomalous Green's 
functions are interrelated. However, here two anoma- 
lous polaron Green's functions are proportional to the 
anomalous one-phonon-cloud propagator Lp{x\x'). Nev- 
ertheless, Dyson's equation for the anomalous polaron 
Green's functions requires the knowledge of the normal 
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X ^ G" x ' X ^ G" x' X ^ G" 
\l G" _ 




1 

"2 Ji 



/> X(/)3(a;,zi,Z2 I ji,j2,x') _ i 



(ii) 



(i) 



(k) 



FIG. 3: More complex diagrams which appear in second order perturbation theory for the normal polaron Green's function. 



polaron Green's function. Then, by summing of corre- 
sponding diagrams we obtain: 

i^pcro^(x, r|x', r') = r2pCTo^(x, r|x', r') 

+ /dri Apcrcr(x, t|xi,ti) 

Xl,X2 

Xt(xi - X2)Fp^7f(x2,'ri|x',T') 

+ X! d.Tinp^^^{-K,T\:>Ci,Ti)t(x2-:><Li) 

Xl,X2 

xGpCT,o^(x',r'|x2,ri) 



with 



l]p(x|x') =F;(x|x')+>S>(a;k'): 



(60) 
(61) 



where Yp{x\x') is the sum of all strongly-connected dia- 
grams for the anomalous Green's functions. This quan- 
tity is analogous to the function Zp{x\x') but differs from 
it by the direction of external electron lines. The Fourier 
representation of Eq. (60) has the form: 

Fp„a:(k\iujn) [1 - e(k) Ap^<^(k|ia;„)] = ripcro^(k|iw„) 

X [1 + e(-k) Gp^(-k| - iLUn)] . (62) 



We obtain by making use of Eq. (59): 

flpaS-mi^n) 



Fpawik\ii^n) 



[1 - £(k) Apcrcr(k|iLJ„)] 
1 

[l-e(-k)Apg^(-k| -iujn)] 



(63) 



The anomalous quantities Clp, flp, Fp and Fp are pro- 
portional to the anomalous one-phonon-cloud propagator 
if, which in the strong-coupling limit is an exponentially 
small quantity. Therefore, it is more important to con- 
sider the propagators Fe{x\x') and Fe{x\x') defined in 
terms of the electron operators a^r (t) and ag. (t) and not 
in terms of the polarons operators c(t) and c^{t). In ad- 
dition we have to discuss the normal electron propagator 
Ge{x\x'). These functions are defined by 

G,ix\x') = -(Taxa(T)ax'a'(^')C^(/3)>§, (64a) 
Feix\x') = -(Taxa(T)ax'a'(T')C/(/3))§. (64b) 

Fe{x\x') = -~{Ta^aiT)a^'a'{T')UiP)yo- (64c) 

Diagrammatical contributions to the first two functions 
are shown in Fig. 5 and 6, respectively. 
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Xips{x,x',ii,i2 \jhj2) Xcp^{x,x', 11,12 I jl,i2) X(/?3(x,x',«i,i2 \jl,j2) 



(e) 



(f) 



Jl. 



,J2 



(g) 

/ \ 

1 / \ 

2 / 













22 ll 








X(/?3(x,x',«i,i2 1 jl,j2) 


X(/?3(x,x',ii,i2 1 ji,i2) 


X(/?(x, x') 


(h) 


(i) 


(k) 



FIG. 4; Diagrams of second order perturbation theory for the anomalous polaron Green's function. The factors tfi„ account for 
the phonon-cloud contributions. 



Figure 5 shows the diagrams contributing to the nor- 
mal one-electron propagator in the presence of acoustical 
phonons (clouds). We again find weakly connected dia- 
grams which can be divided into two parts by cutting 
one electron line like Ci, C3, C5, and Cy . Furthermore, we 
introduce normal, 'E{x\x'), and anomalous, and 
'E.{x\x'), mass operators, of which the simplest contribu- 
tions are shown in Fig. 7. For example, diagram ai is the 
renormalized tunneling matrix element, whereas (b) and 
(c) are the simplest contributions to the anomalous mass 
operators. In analogy to the rectangles representing ir- 
reducible Green's functions, G^"", or non-full cumulants, 

r°, 9^, and Q„ in Figs. 5 and 6, we introduce here the 
correlation functions Zf,(x\x') for the normal state and 
ye(a;|a;') and Ye{x\x') for the superconducting state. For 
example, to (d) and (e) in Fig. 5 corresponding diagrams 
contribute here to Ze{x\x'), while to (d) and (e) in Fig. 
6 corresponding diagrams contribute to the correlation 
function Ye{x\x'), which leads to 



Ae{x\x') ^ G'^(x\x') + ZAx\x'), 

Qeix\x') = F°{x\x') + Ye{x\x'), 
ne{x\x') =F°{x\x') +Ye{x\x'). 



(65a) 
(65b) 

(65c) 



Thus we have introduced the main dynamical quanti- 
ties which determine the to Fig. 5 and 6 corresponding 
diagrammatical structure. This allows us to derive Dyson 
equations for the electron Green's function system. The 
full electron Green's functions can be expressed as 

Ge{x\x') = Ae{x\x') + Ae{x\xi) Ee(a;i|a;2) Ge(a:2|a;') 

- Ae{x\xi)Ee{xi\x2) Fe{x2\x') 

- ^leix\xi)'Ee{x2\xi) Feix2\x') 

- il,e{x\xi) Se(a;i|a;2) Ge(a;2|a;'), (66a) 

Fe{x\x') = fle{x\x') + Qe{x\xi) T.e{x2\xi) Ge{x'\x2) 
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a:' x' X 

(dl) (ds) (d;j) 



Q 



(d4) 



(ds) (de) (e) 

FIG. 5: Diagrammatic contributions to the normal one-electron propagator in the presence of phonon clouds. 



- Qeix\xi)Ee{xi\x2) Fe{X2\x') 

- Ae{x\xi)Ee{xi\x2) FeiX2\x') 

+ Ae(a;|a:i) Se(a;i|a;2) Ge{x'\x2), 



(66b) 



Fe{x\x') = Qe{x\x') + ^e{x\xi) T,e{xi\x2) Ge{x'\x2) 
- fleixlxi) Ee{xi\x2) Fe{x2\x') 

+ Ae(a;|a;i) Se(x2|a;i) Feix2\x') 

+ Ae(a;i|a;) Se(a;i|a;2) Ge(x2|a;') (66c) 

Here x stands for (x, a, r) . Double repeated indices im- 
ply summation over x and a and integration over r. All 
quantities in these equations are renormalized functions 
containing all diagrammatical contributions to the nor- 
mal and anomalous one-electron Green's functions. The 
Fourier representation of the first two equations can be 
written as [k — k, iujn] : 



= GU^k) [1 - A^(-fc) E^(-fc) + nl.ik) S^^(fc)] 
+ F^^ik)[A^{^k)E^,ik) +Tlt,{k) Kik)] , (67a) 

^%Ak) 

= F:^{k)[l ~ K{k) Kik) + KAk) sL(fc)] 

- GU-k) [^%^{k) S^(-fc) + A^(fc) Sl,(fc)] . (67b) 

The solutions of these equations are 
1 



G'A-k) 



dAk) 

-K{k)[K{k)AU-k) + ni^{k)%,{k)\], 

(68a) 

+ S^^(fc) [A^ (fc) A^(-fc) + n%^{k) (fc)] } , 

(68b) 
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F,{x\x') , 
< > = 




(C6) 



(cr) 



i ^ 1 



(C8) 




1 ra^ 



(di) (ds) (d;,) 



X X 

(d4) 




(ds) (de) (e) 

FIG. 6: Diagrammatic contributions to the anomalous one-electron propagator in the presence of phonon clouds. 



where 



%,{k) [AUk) AU-k) + ni^ik) %,{k)] }, 

(68c) 



= i-AUk)Kik)^AU-k)j:u~k) 
+ ni^{k)%„{k) + %^{k)^%{k) 

AUk)AU-k) + nUk)KAk)_ 

K{k)^U-k) + Kw{k)KAk)\- (69) 



relations: 

KAk) - 
Kwik) = 
Kwik) = 

and hence. 



-F§J-k), Ft,{k) = -Fl^i-k), 

-n^^A-k), %Ak) = -KA-k), 

-^A-k), %Ak) = -KA-k), (70) 



dcr(fc) = da-{-k). 



(71) 



The functions F ^ and obey the following symmetry 



These equations for the renormalized electron Green's 
functions are exact. Since they do not contain the ex- 
ponentially small anomalous phonon-cloud propagator 
(f{^W)i superconducting pairing is easier to achieve by 
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Ee{x I x') = 





(c) 



FIG. 7: The simplest diagrams contributing to the normal, Se(a;|a;'), and anomalous, He(a;]a;') and H<;(a:;|a;'), mass operators. 



electrons without phonon clouds but moving in the envi- 
ronment of the clouds belonging to other polarons, than 
by polarons moving in the same environment. We can 
now switch off the superconducting source term, which 
means that F° and F are identically zero. However, 
the functions fi^ and ^ survive in this limit and 
are equal to the order parameters of the superconduct- 
ing state, Y^- and Y-^^, respectively. 



IV. SOLVABLE LIMITS 

The three correlation functions Z^a, Yea.w and Ye^.a 
are the infinite sums of diagrams which contain both par- 
tially and completely irreducible many-particle Green's 
functions. In order to obtain a closed set of equations 
which can be solved (at least numerically), we restrict 
ourselves to a class of rather simple diagrams which, how- 
ever, contain the most important spin, charge and pairing 
correlations. 

One way to do this, is to check how the individual dia- 
grams are influenced by the phonon fields, for example by 
distinguishing the case of moderate coupling, when Eq. 
(26) can be used, from the strong-coupling case where 
Eq. (25) holds. This helps to eliminate from the di- 
agrams the less important ones. For example, in the 
strong-coupling limit (j){x\x') oc (5x,x' holds, which allows 
to discard all renormalized tunneling matrix elements of 
the form t(x'— x) 0(x — x'|0). In this limiting case the 
narrowing of the electronic energy band is maximum, i.e., 
its width is equal to zero. Since this extreme case is a bit 
unrealistic, we will in the following consider the case of 
moderate electron-phonon coupling when also the band 
narrowing is moderate and Eq. (26) must be used. After 
summing the infinite series of the most important contri- 
butions we obtain the result which is shown graphically 



in Fig. 8. 

It is evident that in this approximation for the cor- 
relation functions no closed set of equations is obtained 
because of the complicated nature of mass operators. So 
we have to simplify the latter quantities of which the 
simplest diagrams are depicted in Fig. 7. For the normal 
mass operator we will use the contribution (ai) in Fig. 7 
which is given by 

Eeix\x') ~ i(x - x') (/)(X - x'|0) 

~ i(x-x')exp(-cri(x-x')V2)- 

For simplicity we replace in the exponential function the 
distance |x — x'| by the lattice constant a being a char- 
acteristic length over which the electrons tunnel: 

T,eix\x') ^ i{x - x') S{t - t') 

= i(x - x') exp{-Wp) 6{t - r'), (72a) 
^(Tia-. (72b) 



This result means that tunneling of phonon-fields leads 
to electronic band-narrowing effects by which the bare 
energy e(k) is replaced by e(k) = e{'k)e~^'' . For moder- 
ate electron-phonon interaction the quantity Wp is about 
unity. With respect to the anomalous mass operators, Se 
and Se , we observe that they are smaller than the normal 
one and can therefore safely be neglected. This will be 
used when expanding the equations close to the super- 
conducting transition temperature. 

Another approximation is related to a simplification 
of the exact Dyson Eqs. (68) by omitting all anomalous 
mass operators, which yields [k = (k, ico)] 



■{Al(k) [l~e{--k)AU-k)], 



Dm 

£(-k)r;^(fc)Fl^(-fc)}, 



(73a) 
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FIG. 8: Schematical representation of the renormalization function Zea and the superconducting order parameter Yea,-cr- 
Rectangles with four inner points depict the two-particle completely irreducible Green's functions G2^^ or partial cumulant F^. 
Double lines depict the renormalized one-particle Green's functions of electron and polaron kind, short-dashed lines stand for 
the hopping matrix elements and long-dashed lines represent the direct electron-electron interaction V {i — j) . 



(73b) 



X ^ ^ / dTidT2G'2*''[CT,T;CTi,ri|cr2,r2;cr',T'] 

X t(j -x)i(x-i)0(O|ri - T2) Gp(i, (72, r2|j, (Ji, n), 

(74) 



Dm 

D%{k) = [1 - £(k) A^(fc)] [1 - e(-k) A^(~fc)] , 

+ e(k) £(-k) Y^m%A-k), (73c) 
A^(fc) = G°(fc) + Ze,(fc). (73d) 

These equations are identical in form with the Dyson 
equations for polaron superconductivity mediated by op- 
tical phonons^^. The difference to the previous work is 
related to the appearance of the renormahzed energy e(k) 
and new correlation functions shown in Fig. 8. These 
irreducible functions depicted by rectangles are on-site 
quantities with equal site indices. Hence, all right-hand 
parts in Fig. 8 are proportional to 5x,x' meaning that 
Zeaa and l^eo-CT are also local functions and correspond- 
ing Fourier representations to be independent of the po- 
laron momentum k. In the diagrams x and ii stand for 
(x, cr, r), and (ii,cri.Ti), respectively, whereby summa- 
tion over ii, ji, 12, J2 and cri, (T2 and integration over ti 
and T2 is assumed. These quantities have the following 
analytical structure: 

Ze(x, a, t|x', a', r') - -i^x.x' ^'(^ " 

j 

X / dTidT2r3((T, r;CT',r'|ri,T2)»^^''-'5x,x' 
Jo 

V V / dTidT2G'^"[CT,T;CTi,ri|cr2,r2;CT',T'] 



(Tl ,(T2 I] 

X i{j -X.) f(x-i)Ge(i,cr2,r2|j, (Ti,ri) - (5x 



y^(x,a,T|x',-a',r') = -i<^x,x' ^^i^^2 

o" 1 ,f2 ii ,12 

X G2 ''■[cr, r; -fj', rVi, n; 0-2, r2] i(x - ii) i(x - 22) 

X F''(n,0-i,ri|z2,cr2,T2) - ^f^x.x' 

X J2 Xl'^2*''[cr,-^;-CT',^Vi,n;CT2,T2]t(x-n) 

CTi,(T2 1112 

X f(x- i2)-F"'(«l,O'i,Ti|i2,Cr2,T2)0(O|Ti - T2) 

X 0(O|r2 - n), (75) 

and corresponding Eq. for Y^a-.a- Since Eqs. (74) and 
(75) are the result of summing an infinite series of dia- 
grams, the thin lined representing one-particle propaga- 
tors are replaced by full normal electron (Ge) and po- 
laron (Gp) and anomalous (F^) functions. Fourier trans- 
formation of these quantities leads in case of spin-singlet 
channel of superconductivity to 

Z,.(zc.) = ~iy2n'^r..(^c.)-^ ^ [e{k)r 



X Geai ,ai (k| icji ) G" " [ct, iuj;ai, iuji \ ci , iui ; cr, iu] 
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X G2^'^[a,iuj; (71, iuJi\cr I, iuJi; a, iu>], (76) 

Yea,-a{ii^) = -^^ e(k) e(-k) Fg^i ,-<7i (k|ia;i) 
X G2^^[a, iio; —a, —iu)\aijiu)i; — cri, —icui] 

(Ti,k a;ir2iri2 

X Fea,,-ai (k|i(a;i_Qi + ^2)) (t>W 

X G2 — iwlcijiwi; — (Ti, — iwi], (77) 

where 

k 

/■/3 

r<x,a'(T-T') = / rfridT2r3((T,T;f7V|ri,T2), (78b) 
JO 

G2 [(Ti , iti^i ; (72 , 1 0-3 , ; (T4 , «W4] 
= /3 (5(aji + 0^2 — W3 — CJ4) 

X G2"'[a-i,ia;i;cT2,«W2|a-3,«W3;cT4,«(wi + W2 - W3)]. 

(78c) 

I 



where is the Kronecker symbol for Matsubara fre- 
quencies and 

X{iuj) = + /z, A(iw) = iio + ji — U, a = —a. (81) 

For the present study and [/ in Eqs. (80a-c) and 
(81) are the renormalized quantities of Eq. (13). The 



From Eq. (76) for Ze{iio) we obtain the following expres- 
sion for Ae(zw), 

UJl ,k,(Tl 

X Geo- (k| iuJi ) G2 [cT, zo;; (71 , icJi | ai , icJi ; tj, ia;] 
- ^ E ^ E (k) Gp.(k|i(a;i + Oi)) </)(ifii) 

k,(Ji oJi^i 

X G2''^[a,iuj;ai,iu!i\ai,iLiJi;a,iuj]. (79) 

Here, the renormalized and unrenormalized tunneUng 
matrix elements accompany the electron and polaron 
propagation, respectively. Equations (77) and (79) to- 
gether with the expressions for the onc-particlc Green's 
functions and the definitions of the irreducible Green's 
functions in Eq. (50) and Kubo cumulants in (55) de- 
termine completely the properties of the superconduct- 
ing phase and allow to discuss the influence of strong 
electron-phonon interaction. 

The three irreducible two-particle Green's functions al- 
ready calculated (see Eqs. (5-7) in Ref. 13) are given by 



(80a) 



(80b) 



(80c) 



fore-standing equations are generalized Eliashberg equa- 
tions of strong-coupling superconductivity for the case 
that strong electron correlations have been taken into ac- 
count in a self-consistent way. In spite of the approxima- 
tions involved the equations are rather complicated. In 
order to gain further insight into the physics behind Eqs. 
(77) and (79), we will linearize the equations in terms of 
the order parameter l^,^, but not in terms of A^. Then 



G2'^^[cr,itJ\<J,ioJi\a,iLOi\<j,iuj\ = 



PU^[1 - ,5^,^,](1 + e^^)(e^^ + e^(^^-^)) 
Zq Xiiuj) Xiiu) X{iu!i) X{iLOi) 

U( pUe^Pf^ie-P^ ^l) 



G2[a,itJ\G,itJ\\a,iw\\a,ito\ = — . _ _ _ _ 

Z[) y Zi) X[iuj) X{iuj) X{iuji) X{iu)i) Xyiu)) Xyiuj) X{iLOi) X{iu}i) 

{e^i^^^-u)-l)\X{^u) + X{^o^l)\ , 1 



+ 



[A(ia;) -|- A(za'i)] A {iuj) X {iuii) 

1 1 

+ 



1 1 
+ 



A(iw)A(iw) \X'^{iuJi) A^(iwi) 



X{iu)i) X{iuJi) \'^^(^'^) X^{iw)J X{iu)) X{iu)) X{iu)i) X{iu)i) 



G2 "^[cr, iw; C7, — ia;|cr, iwi; (T, — iwi] = 



I3U 



Z<,\[n^-{iu:mn-Uf-{iu:f] 



^0 



-/3f/\ 



U 



(l + e^") 



+ 1 



u 



2m - Uj [(m - C/)2 - (iu;)2][(M - [7)2 - (ico,)^] 

I 
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the critical temperature Tc of the superconducting tran- sition can be obtained from 



1 e(k)£(— k)y^CT,CT(*wi)G2'''[o')*w;(T, — ia;|c7, iwi;c7, — iwi 



k,aJi 

—y-y 



[1 - £(k) - £(-k) Ae^(-ia;i)] 

e(k) £(-k) yeo-,?f(«('^i - ^^1 + ^2)) H'i'^i) 4'i'i^2) G^^'^la, ioj;a, -iuj\a, ioJi;a, -iuii] 
[1 - £(k) Asaiiooi - ifil + iO.2)] [1 - e(-k) Aea{-ioJi + if^i - iO.2)] 

£(k)]2 Ae<^i (iwi ) [a, iuj;ai, zwi | tri , iui ; cr, iw] 



1 - e(k) Ae<., (iw) 



E 



G2 [(T,ico;ai,tiO-i\ai,icoi;a,iui\. 



k,aJi,ni,(Ti 



1 - £(k) Ap„^ {iu) 



(82) 



(83) 



In order to determine Tc it is necessary to solve Eq. (83) 
for Aea{iui) and to insert it into Eq. (82). The next our 
approximation is the omitting in the right-hand part of 
eq.(83) of the polaron fimction Ap cr.The reason for such 
approximation is the presence in this term of the phonon 
-cloud propagator (^(ifii) which makes the appearance in 
the denominator of large quantity Wc- Since coc is larger 
than other typical energies involved, the term under dis- 
cussion is at moderate coupling strength smaller than all 
other terms in Eq. (83) and can therefore be omitted, 
which leaves 



1 

1 

w 



1 [e(k)]^Ae(iwi) G2*'"[ct, iw; a, iwi |(7, iwi; cr, iw] 



1 -£(k) Aea(«w) 



\p |£(k)pA^(ia;i ) G^ [a, iu);a,iLOi\a,iu}i; a, ito] 



k,u;i 



1 - £(k) Aea(iw) 



where 



(84) 



(85) 



Equation (84) is identical with the corresponding 
equation for the single band Hubbard model without 

phonons^'' if we replace in Eq. (84) the renormalized 
quantities /x, U, e(k) and G^ {iu>) by the initial quan- 
tities /io, Uq and £(k). This means that we have re- 
duced the investigation of superconductivity in frame of 
the Hubbard-Holstein model to the analogical problem 
with respect to the single band Hubbard model. 

It is instructive to analyze the contributions from the 
two spin channels by considering the quantities 

(■ \ 1 [£(k)]2Ae.(lC^l) 



£(k) AeaiiiOl) 



X G2 [a, icj; a, iuJi \a, iuJi ; a, iuf\ 

Jn 



. ^ 1 Y-v- [e(k)]2Ae5^(ia;i) 

uJi k 



£(k) AeCT(iWl) 

X G%^'^[a,iw-,a,iwi\a,iuJi-,a,iw], 
and using the notation 



' N ^l-s{k)A,,{iuj) 

= 1 V - '^^^ 
iv 4^ 1 



£(k) AeaiiLo) 



(86a) 



(86b) 



(87) 



Here it is assumed that e(k) = e(— k) holds with 
^,^£(k) = 0. We replace sums by integrals, 



(88) 



ttW 



- (2e/Wy 



1, |£| < 1^/2 

0, |£| > W/2 



(89) 



W is the renormalized band width W — We^^p and pq 
the semielliptic model density of states. Since we do not 
consider magnetic solutions, the spin index can be omit- 
ted. Making use of (80a-b) for the irreducible functions 
we obtain 



X+(«w) = 



X{iLu)X{ioj) 



(3X{iuj) X{iuj) 



(90a) 



X-(ia;) = — 
^0 



U ( 2(/x - [/)(e'5(2A'-^) - 1) <l>%{-iuj) 



P{2ij.-U)[X{iLo)X{iLo)Y 
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-U 



+ 



[X{iuj) X{ioj)f 



(l + e/3M)(0^ + 20) 



1 



where 



1 



P^X{iu;^)X{iw^y 

-0^'^'^^'^^^{x^)-Jik) 



(90b) 



(91a) 



U 



/3'5[A(^'^i)A(^c.i)F' 
,. _ 1 Xjiuji) + Xjiuj) -tJj^+^^^o(l)l{iu;i) 



(91b) 
(91c) 



and V'i^+aji,o = 1 — <5w+a)i,o- In case of half fiUing when 
/i = 17/2, X{iuj) = ioj + fi, A(iw) = ioj — fi, wc have the an- 
tisymmetry property, (jf{—i<jj) = —(jf{iuj), and therefore 
(j) = =0. Furthermore, we find at half filling. 



X+(iw) = - 



[(iw)2-/x2]2' 
2M2(/)'=(iw) 



= 2x+(«w), 



(92a) 
(92b) 



and therefore Eq. (84) is equal to 

X+(iw) + x-(«w) = 3x+(ia;). 



(93) 



Away from half filling we will for simplicity omit the 
contribution of the anti-parallel spin channel in Eq. (84) 
and introduce instead a correction factor /g, which is 
three at half filling and diflferent from three at general 
filling. The final equation to be investigated is then 

/, ^ [£(k)]2Ae(ic.i) 



pN^ l-e(k)Ae^(ia;) 

X G2^'^[cr,iuj;a,iu!i\cr,iuJi;cr,iLo]. (94) 

Further simplifications can be made regarding Eq. (82) 
for the critical temperature. We note that the main con- 
tribution to the last term results from the minimum val- 
ues of frequency difference f2i — Cl2- When Oi = 0.2 we 
get after summing over fli: 



l5^[<A(ifii)P = -^cothi/3a;e + 



so that 



2sinh2 ^jSuj, 



; (95) 



f c = 1 + COth i/^Wc H 



(96) 



can be used as a common factor in the remaining function 
for the superconducting order parameter, which is 



f{iuji) G2 *'"[cr, iui; fj, - 



k,aJi 



[1 - £(k) Kea{i^l)][l - £(-k) Ae^(-lWi) 



(97) 



In order to solve Eq. (97) for Tc, we introduce a new 
function, 



e(k)£(-k) 



N^[l- f(k) Ae.(iw)][l - £(-k) Ae^(-i0 



(98) 



Aeo-(ia;) - Aea(-iw)' 

which allows to rewrite Eq. (97) in the form 

f 

X G2^'^[a,iu};a,—iu:\a,iuJi;a,—iuj\]. (99) 
Using furthermore 



_ 1^ ^ (l)'"'{iL0)Yea;a{iuj) 



(100a) 



^ 1^ 0*'=(iw) i;<T,o^(«w) 



c = 



^Z^(^_t/)2_(,^)2' 



Q(ia;) = 1 + 

3WS to 
^e<7,a(«w) 



^0 



allows to obtain the solution: 

C//e[l + C//(2/x-C/)](l + e'3'^)Vi 



(100b) 

(100c) 
, (lOOd) 



Zo Qiiuj) - (tu;)^] 

Ufc [1 - U/{2(i - Z7)](e^^ + e^(2^-^)) V-a 
Zo Q{iij) [{fi - [7)2 - (ia;)2] 

(101) 

where symmetry properties of Yea.aii^) and 0*'^(ia;), 
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has been used. The two constants tpi and V-'2 are deter- 
mined from the foUowing system of equations, 



1 + 



'ip2 



Vll fc 




J" 7" 

U 




^0 


^/i — 




'712/c 


[l- 


J/ 




^0 


2/i- 


u 


'722 /c 








^0 




u 


»7l2/c 




U 




^0 




u 



Where rjij are given by the summations of 



= -y 



ryi2 = 



P ^ Qiioj) - (jw)2]2 ' 



a ^ Qiioj) - - [7)2 _ (i^)2] ' 



U 



'?22 = 7^E 



/3^QM [(/,-[/)2_(ia;)2] 



212 ■ 



(104) 



The critical temperature is then obtained by setting 

T = Tc in the equation: 



1 + ^11 + 

^0 



U 



1 + ^11 



1 



X (1 + e'^") = 0. 



(l + e'^") 

(105) 



U 



Equation (105) is invariant under the particle-hole trans- 
formation, 

ria ^ 1 - ricr, 

[/^ -(2/i- [/), 
Zo(/x) ^ Zo(/x)exp[/3(;7 - 2^)]. (106) 



V. METALLIC, INSULATING AND 
SUPERCONDUCTING PHASES 



In order to obtain a better understanding of the prop- 
erties of the main equations, we will first check whether 
the normal state determined by Eq. (94) is metallic or di- 
electric. This can be done by analyzing the renormalized 
density of states, p{E), which is given by 



1 



p{E) = luig{E + iQ+), 



5W = ]vZ.T^ 



A(ia;„) 



£(k)A(ia;„)' 



(107a) 
(107b) 



where A(iw„) has to be calculated from Eq. (94). The in- 
tegration over k is again done by using the semielliptical 
form of model density of states in Eq. (89). This gives 
for the quantities (87) and (107b) the following result^'': 

(p^iiuj) = ^W<fi%ioj), (108a) 
ip%iij) = (1 - [1 - A.^{iuj)]^/^f/A.^{iij), (108b) 

g{i(j) = {A/W) (1 - [1 - A2(ia;)]V2)/A(j^), (lOQa) 
l{iuj) = \WK{iuj), (109b) 

yielding for the renormalized density of states (DOS) in 
Eq. (107a) 



ME) 
wW 



Im — 

wW 



1-A2(£;) 



HE) 



(110) 



where A{E) is the analytical continuation of A(iLUn). 

We now address Eq. (94) which determines the normal 
state of the system. By Using (90a) and (108) we rewrite 
it in the form: 



Aiiu) = iVKG^(^^) - lmx+{tLo), 



(111) 



with x+{i'^) given by Eq. (90a). In the limit E —^ Eq. 
(Ill) becomes 



A 



a 

A4 



-■\/l-A2^ 



(112) 



where A = A{E + iS) for E = 0. The two parameters a 
and b and the function (0) are given by 



C/2(l 



o/3M)(e/3M 



p/3(2p-l/) 



b=WG''{Q)-\fs 



|3WU^{l + eP^'){eP^' + 



„/3(2m-C/) 



G^(0) 



1 



+ 



m(m - u)zl 



p. p. — u 



n ■■ 



1 



1 + e^^ 



(1-n), 



(113a) 
(113b) 

(113c) 
(113d) 



with d) given by Eq. (91). The first term in Eq. (113c) is 
the value of the local one-particle Green's function of the 
Hubbard model for £^ = while the second term is the 
corresponding contribution from the intersite Coulomb 
interaction discussed in Appendix B, see (B16). 

Equation (112) has been discussed for the simpler case 
of half filling when 6 = in Refs. 17 and 27. Away from 
half filling we have n ^ 1, fj, U/2 and 6^0. As 
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Eq. (113a-b) show, the parameters a and h depend on 
chemical potential ^, strong-coupling parameter U and 
mean electron number per lattice site n. Moreover, b 
depends on the intersite Coulomb repulsion. Although 
parameters a and b are not independent of each other, 
we first try to get information in case that 6 = 0, i.e., for 
^ = G^{v) = 0, which holds for K{-iu)) = -k{iuj). In 
the half-filled band case when ^ — U /2 and n = 1, the 
quantity a is equal to 



(114) 



The equality b=0 allows to determine A in Eq. (112) from 
a simpler relation, 



1- Vl-A2 
A2 



For a > 1 the first factor gives A^ 



= 0. 



1- Vl- A2 
A2 

''(115) 
= a(2 — a] and hence, 



A = ±^a(2-a), a < 2, (116a) 
1 = ±i^/Ma~~^, a>2. (116b) 

By inserting these solutions in Eq. (109) for the renor- 
malized DOS we obtain results different from zero only 
for the expression with lower sign in Eq. (116b), i.e., for 
a > 2. Hence we obtain a metallic state at half filling only 
if the Coulomb interaction is less than a critical value* , 



(117) 



In this case there is no gap at the Fermi level and the 
renormalized DOS becomes (a > 2) 

p{0) = A/iirW) r(0), r(0) = ^/ {a - 2)/a. (118) 

The insulating (dielectric) phase exists if the inverse con- 
dition holds, a < 2, U > Uc, leading to the opening of an 
energy gap at the Fermi level. 

Away from half filling, for a > 2 and b sufficiently 
small, the solution of Eq. (112) can be obtained from a 
series expansion in powers of 6, 



Aft = A 



(a -1)6 (a2-5a + 3)62 



2{a-2) 8A(a- l)(a- 2)2 



(119) 



where A is given by Eq. (116) in zero-order approxima- 
tion. This leads to 



r(0) = 



a - 2 {2a-^ - Ba^ + 12a - 5)6^ 
^ 8a3/2(a_ i)2(a_ 2)5/2 



(120) 



This result shows that not b but b/{a ~ 2) should be 
taken as expansion parameter. Therefore the expansion 
is not correct very close to a = 2. Another peculiarity of 
Eq. (120) is its even character in b, which follows from 
the fact that in the solution of Eq. (112), A, changes 
sign when the sign of b is changed. Therefore, if we have 



the solutions Ai±iA2 of Eq. (112), then correspondingly 
— Ai ^ iA2 are solutions for 6 < 0. The different signs of 
the real parts do not matter since Eq. (110) depends only 
on the absolute value of Ai . The numerical investigation 
shows that for a > 2 the role of b is not decisive, i.e., 
the system remains metallic. However, in the range < 
a < 2, for which the system is insulating when b — 0, the 
influence of b is decisive because there exists a critical 
value, bc{a) ~ 1 — a/2, such that the system is metalHc 
for |6| > be and insulating for |5| < be- Note that the 
parameters a and b are not independent each other and 
therefore such extremal parameter values as a small and 
b large or vice versa are not admitted by their definition 
(eq.113). 

The metallic state exists for low and high band filling 
(small and large values of fi) and near half filling, (/i ~ 
U/2), provided that U < Uc- The physical role of b is 
to enhance in each case the tendency towards metalicity. 
There is also the case of reappearance of the metallic 
phase for a < 2 and |5| > bc{a). The variation of bc{a) is 
shown in Fig. 9. 




Parameter a 

FIG. 9; Variation of the critical parmaeter bc{a) in the range 
< a < 2. For l&j > be is the system metalhc while it is 
insulating for |fe| < be- For a > 2 the system is metallic 
regardless of the value of b. 

We discuss now the superconducting phase transition 
with Tc obtained from Eq. (105). In order to gain some 
insight we need to simplify (j)'"^(iu!) from Eq. (98). In the 
asymptotic limit |ct;| — > cxd this quantity is equal to 



0-(,^) ^ i(M//2)2 



(121) 



whereas in the low-energy limit £^ — > required here, we 
obtain from Eqs. (90a) and (111) the value 



(t^'^iE + id) = (1/a)^ (IF/2)2 



(122) 



The two limits coincide for a — 2. Really we are inter- 
ested mainly in the low energy limit (eq.l22). Because 
of the fast convergence of the sums in Eq. (104) which 
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contain (t>'"^{iuj) and which determine the parameters rjij 
of Eq. (105), wc can rewrite 77^ as 



where 



(d2 - z^) 



??22 



??12 = 



/3 

Uct)"" sr^ 1 
^ o77 



/3 



<3i(^) 



l! 



r _ 1 V 1 



1 



/3^gi(z),(z2-c2)^ 
1 



(123a) 

(123b) 

(123c) 

(124a) 
(124b) 
(124c) 
(124d) 



and 



/,,(e'^" + l)(e'3A' + e/3(2M-t^)) 



2 2 



(125) 



In case that U > the quantity C</''"^ = 7^ is positive. 
In the other case of very strong elcctron-phonon interac- 
tion we have U < 0. In this case we will use the notation 
^0*^= = —0^ (with negative (). Although U < seems 
not to be the generic case, wc believe that in the metallic 
phase in case of strong electron-phonon interaction an ef- 
fective attractive interaction is possible, because the local 
and intersite Coulomb interaction can be largely screened 
by the electron-ion interaction. However, calculation of 
self-consistent screening in not the aim of this paper. We 
just require charge neutrality and take C/ as a parameter 
which allows to discuss the case of effective repulsive as 
well as attractive interactions. We first discuss the case 
[/ > and then U <0. 

The sums in 7„ have been evaluated by contour inte- 
gration with the help of Poisson's formula, see Eq. (C2). 
For the case 7"^ > (c^ — d^)^/4 we obtain the following 
equation which determines Tc- 



0=1 + 



^^4 + I ^474 - (c2 - (Pf 



1 + 



U{1 



Zo{2n-U) 



0.2 sinh j3a\ — oi\ sin (5ai 
cosh /3q!i + cos /Ja2 



+ 



J7 + 



(2/x-;7)(l-e'5(2/'-c^)) 



U\a\ sinh/3ai + sin/3Q!2] 



cosh /3q!i + cos /3a2 



+ 



+ 



tanh i/3c? 



(/x-C/)(e''^ + e''(2^-^))- 



/utanh 5/3c 



(1 



2 t/V(M-^)(^</'"")'(l + e^^)(e^^ + e^^^^~^^) / ^-j^ - {c^ - d^ f 1 

78^2 + c2d2 \ ^474 - (c2 - d2)2 (c2 - ^2) 

/tanhi/3d tanh i/Jd^ a2 sinh /?qi - a-i sin /3a2 sinh^ /3ai + sin^ /3a2 



d 



cosh /3ai + cos /3q;2 



[cosh/fei + cos/?a2]^ 



-— tanh tanh / tanh tanh a, sinh/3ai + a2 sin/3a2 

- V7^ + c2d2 ^ — + — + ^ — — 

d c \ d c / cosh fjai + cos pa2 



where the parameters ai and a2 are given by Eq. (C3). 
For the special case of half filling when /i = {7/2 > and 



r 



^ -^^ 6/3" + 1 ' 

Eq. (126) is of simpler form. 



(126) 



(127) 



= 1 + 



3^4^^4_^^4 1 ;x2 



/3/Lt \ q;2 sinh /3ai — ai sin /3q;2 _|_ ^ "^i snih. j3ai + q;2 sin /3a2 



e/3A» + 1 y cosh /3ai + cos /3a2 



cosh /3ai + cos (3a2 
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- Vt^TT? tanh 



4/,^^° / sinh^ 0ai + sin^ /jaa a/^EZ 



tanh^ 



7^ / d tanh a2 sinh/3ai —«! sin /3q2 

2/z \d/Lt /X y cosh /3q;i + cos /3q;2 



2 , -I - ai sinh /3q;i + 02 sin 

— tanh — 

II cosh (jai + cos pq;2 



(128) 



where 



ai =-^(x/7"±/z4 + ;x2y^', (129) 

If we consider at half filling in addition P/j, ^ 1 and 
Pojc ^ 1 we get 7^ = /x^ and fc = ^ and instead of Eq. 
(129) we obtain 



and Eq. (128) for T = Tc becomes 



(130) 



+ l + 2^\/2+l ^l + 2\/2 + 2^\/2 + l^ =0, 

(131) 



which can not be fulfilled. This shows that for the half- 
filled band case and positive renormalized Coulomb in- 
teraction s-wave superconductivity is not favored. This 
may change for d-wave superconductivity, which would 
be a technically more demanding investigation in view of 
our expansion around the atomic limit, because of the 
explicit k dependence of the order parameter, which has 
to be retained in case of d-wave superconductivity. 



In the following we will discuss the opposite case, i.e., 
C*/**^ < (c^ — d^)^/4, which includes the possibility to 
consider the negative values of ( and hence (^(j)^'^ = —9'^ < 
0. This case corresponds to f7 < 0. The corresponding 
equation for Tc is then 



0=1 



46*4 



^(C2 _ d2)2 + 45(4 



J7(r'**-''"''' - l) \ ft&nh.\j3zi tanhi/3z2 
Zo{2n-U) 



Z2 



+ u[ u + 



(2/i- J7)(l - e^(2^-^))\ /tanhi/3.21 tanhi/302 



^0 



+ 



Z2 



+ W[{ii-U) 



+ 



tanh \(3d {e^^^ + e/3(2M-!7)) tanh 12/3c (1 



-II- 



Zn 



-2 



c 
1 



+ 



Zl 



+ 



/2t/V(At-C^)(</''')^(l + e'^^)(e'^^ + e'^(2M-c/)) l^tanhi/Jc ^ tanh i/3dN^ ^ tanh ^ tanh i/3z2 

/ tanh h(3zi tanh h(3z2 tanh i/3c tanh ^f3d 



Zl 



(c2 - d2) 



\ z\ Zl c d 

^tanh i/Jc tanh i/3d^ ^tanh ^pzi tanh i/3z2 ^ (c^ - (/2)2 + 261 



Zl 



Z2 J ^(c2 - d2)2 + 45(4 



(132) 



where In case of very small 6 values corresponding to the 

^ condition 

Zl =^ (c^ + ± ^(c2 - d2)2 + 4^4 y ' , (133) 

= \/2V ^ ; ^ ^ ln/3|M| =ln(/3i|C/|) >ln3 (135) 

For half fiUing, U < and 6^ < c^d^t wc have Eq. (132) simphfies to (6* ^ 0) , 



1 — '\p-P\U\/'^ ( / o 

C<^- = -e' = (134) 1 - /,^^«^ (i - ] + 

for sufficiently large |[/|, i.e., In |[//3| > ln9. - ^fl IJ'^ fiMA^ - Aj) = 0, (136) 
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where 



1 d" /tanhi/3;u 



(137) 



Because is assumed to be larger than 3, the coef- 
ficients in Eq. (137) can be approximated by 

^1 - ^2 ^ A3 ~ (138) 

2|/z|3 8|^|5 16|/[i|7 



allowing to replace Eq. (136) by 



(139) 



For negative chemical potential, /j. < 0, and /c = 1 + 
l/(/3wc) the equation for the critical temperature has the 
simple form: 



- 46i^ + (2% - 191)i + 24y = 0, (140) 



with 



yielding for small t the solution, 

y 46 t 



t ■ 



where 



1-y 191(l-y)3' 
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^=191^' 



(141) 



(142) 



(143) 



with the requirement that y <C 1 is fulfilled. In the sim- 
plest approximation Tc is of the order. 
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(144) 



showing that Tc is proportional to the rcnormalizcd 
chemical potential in Eq. (14) and increases linearly with 
increasing strength of the electron-phonon coupling pa- 
rameter. 



VI. CONCLUSIONS 

The interaction of correlated electrons and acousti- 
cal phonons has been discussed by using the canonical 
transformation of Lang-Firsov which results in mobile 
polarons consisting of electrons surrounded by the acous- 
tical phonon fields (clouds). A kind of generalized Wick's 
theorem is used to handle the strong Coulomb repulsion 
between the electrons emerged into the see of phonon 
fields. 

In the strong-coupling limit of the electron-phonon in- 
teraction chronological thermodynamic averages of prod- 
ucts of acoustical phonon-field operators are expressed 



by averages of one-cloud operators. For the normal one- 
cloud propagator the Lorentzian form in Eq. (25) while 
for anomalous one the Gaussian form in Eq. (28) has been 
found. Because the latter propagator is considerably 
smaller than the first one, we find that the anomalous 
electronic Green's functions are more important than the 
corresponding polaronic functions. So the superconduct- 
ing phase transition is determined as usual by the appear- 
ance of electronic Cooper pairs, i.e., the pairing of elec- 
trons without phonons-clouds is easier to achieve than 
the paring of polarons with such clouds. 

For the system of renormalized electronic Green's fimc- 
tions in Eq. (64) the diagrammatical structure is ana- 
lyzed and the Dyson equations have been derived, see 
Eqs. (66-69). Besides the full Green's functions the equa- 
tions contain also three correlation functions and three 
mass operators. These quantities have been calculated 
by summing infinite series of diagrams after performing 
appropriate approximations of Eq. (68). Resulting Eqs. 
(77) and (79) for the superconducting state are then lin- 
earized in terms of the order parameter Y^^^-^ leading to 
the final Eq. (105) which determines the superconducting 
transition temperature Tc, which is invariant with respect 
to particle-hole transformation. Further analysis shows 
that the problem of superconductivity in the frame of the 
Hubbard-Holstein model is analogous to the discussion of 
superconductivity in the frame of the single band Hub- 
bard model with appropriately renormalized paramaters. 

For further discussion the normal state properties 
given by Eq. (94) are investigated with respect to the 
metal- insulator transition. For half filling, i.e., one elec- 
tron per lattice site and iJ = C//2, the model yields a 
metallic state provided the renormalized value of U is 
smaller than \^/2W, where W is the electron energy 
band width, which is narrowed by the effect of the phonon 
fields, see Eqs. (112) and (115). The second parameter 
b in Eq. (112) determines the deviation from half- filling. 
It has been shown that away from half filling a metal- 
lic state is favored for U > \/3/2W and b larger than a 
critical value be- 

The search for superconductivity has then been per- 
formed on the basis of Eq. (105) for s-wave supercon- 
ductivity for two cases. In the first case, at half filling, 
(/i = U/2) with positive U, Eq. (105) has been reduced 
to Eq. (131), which has no solution. In the second case of 
negative U the non- trivial solution in Eq. (144) has been 
found. As mentioned before, in a real solid with corre- 
lated electrons the quantity U should be replaced by an 
effective screened parameter. For an overall attractive in- 
teraction mediated by the phonons the Hubbard-Holstein 
model can have a superconducting solution although the 
bare value Uq can still be substantially large. 
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APPENDIX A: LAPLACE APPROXIMATION 

In this section we provide calculational details of the 
Fourier representations ^(if2) and ^(ifi) defined in Eqs. 
(24a) and (24b) by making use of the Laplace method of 
approximation^^. In the strong-coupling limit the inte- 
grand is maximal at the end points r = and r = /? and 
is considerably smaller at other points of the interval of 
integration (0,/3). Therefore, we can replace the initial 
integral in Eq. (24a) by one in which the region of inte- 
gration is limited to the two small intervals (0,to) and 
(/3 — To, /?). In these intervals insertion of the expansions 
(19) lead to 

c0 



Jo JI3-T0 



[ dre*"^-'""=('^-^). (Al) 



Then, because in the strong-coupling limit the collective 
frequency Wc is large, we can replace tq by infinity, which 
yields a Lorentzian for (Al), 



(iO)2-(w,)2- 



(A2) 



If we take into account the space dependence and expand 
in terms of small distances, |x|, we obtain: 



C7(x|r) « a(0|T) - iaix2, (A3) 



where 



1 

6N 



coth ^/3wk 



sinh ^/?Wk 



(A4) 



For small values of |x| the function in Eq. (21) can then 
be written in factorized form. 



</.(x|r) </.(x)0(O|t), 



where 



<^(0|t) = e-(o|o)+-(o|-) = 1 ^0(jn)e 



-iQr 



</.(x)=e---V2 = ^^^(g)e-iqx 
n 

^(9) = (27r/ai)3/2e-«'/-^ 



(A5) 



(A6) 



(A7) 



(A8) 



For the Fourier representation of the anomalous 

phonon-cloud propagator (/5(x|t) we consider first |x| = 0. 
In this case we need the r expansion of a (0|t) near the 



midpoint of the interval t = (5/2 where this function is 
minimal. Near this minimum we use the following ex- 
pression: 

a(0|r) a{0\(3) + l<j"{0\f3/2) (r - (3/2)\ 

This approximation can be used in the integral (24b), 
which allows to rewrite it as 



P/2 + T0 



^.^g-cr(0|0)-a(0|/3/2) 



13/2-To 



X ginr-^(r-/3/2)V'(0|/3/2)_ 



(A9) 



The width of the small interval, 2to, is now extended to 

infinity because the second derivative, (t"(0|/3/), is large 
in the strong-coupling limit, which yields a Gaussian dis- 
tribution. 



ip{i^ln) = \/27r/cT2 e 
(72 = (t" (01/3/2) 



-<7(0|0)-cr(0|/3/2)-Fi/3a/2-a2/(2(72) 

(AlO) 



This allows to get an approximate expression for Eq. (38). 
With Eqs. (30) and (32) we have, for example, to evaluate 

(f>2{xi,ifli;x2,ifl2\x3,ifls;x4,ifl4) 



X exp {a(xi - xalln - ra]) + (t(xi - X4||ri - T4I) 

-I- a(x2 - X3||r2 - T3I) + C7(x2 - X4IIT2 - Til) 
- Cr(xi - X2IIT1 - T2I) - C7(X3 - X4IIT3 - Ti\) 

-2cr(0|0)}. (All) 

Equation (All) leads to a sum of 24 fourfold integrals 
with different chronological order of t„, n = 1 . . .4 (for 
example, ^Ti>r2>r3>r4 jg defined by from /3 > n > T2 > 
T3 > r4 > 0), of which only 16 make an essential contri- 
bution in the strong-coupling limit. It is convenient to 
combine the 16 terms pairwise like 

[n > T3 > T2> T4 & T2 > T4 > Tl > T3] , 

[n > Ts > T4 > T2 & T4 > T2 > Tl > T3] , 

[n > T4 > T2 > T3 & T2 > T3 > Tl > T4] , 

[ti > T4 > T3 > T2 & T3 > T2 > Tl > T4] . 

Further pairwise terms are obtained from the chronologi- 
cal orders by changing in the first two groups the order of 
Tl and T3, and in the last two groups the order of n and 
T4,. All integrals arc then calculated by using the maxi- 
mum possible value of S in Eq. (32) in the Laplace ap- 
proximation, i.e., maximal contributions arise from pos- 
itive terms in Eq. (32) and coinciding arguments in each 
a function. In addition the t space, for which the inte- 
grand is maximal, should be large enough. For example, 
in the integration over the first pairwise terms we take 
|ti — T3I and |t2 — T4I as well as |xi — X3I and |x2 — X4I 
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as small quantities. This leads with ti — = ti and 
T2 — T4 = t2 and, assuming for simplicity, Xi = X3 and 
X2 = X4, to the following argument of the exponential 
function in Eq. (All): 

a{0\ti) + (T{0\t2) 

+ (t(xi - xaln -T2 + 12) + cr(x2 - xi|ri - T2 - ti) 

- (j(xi - X2|ri - T2) - ct(xi - X2IT1 -T2-ti+ t2) 

-2(j(0|0), (A12) 

which simply reduces to —uic{ti +12) when expanding in 
ti and t2- The corresponding integrations over t\ and t2 
in the interval (0, tq) — > (0, 00) lead to 

^ri>T3>T2>T4 

= / dri / (^7-2 e''^i^^i~^^'+''^^^^^~^*^ 
Jo Jo 

Jo Jo 



^Xl ,X3 ^X2 ,X4 



(-Wc + iri3)(-Wc + 1^4) 

< fdTi I d7-2e'"("i-"^)+'^^("^-"*). 
Jo Jo 



(A13) 



^r2>T4>ri>T3 ^jgg^g f^Q^ ^ri>r3>r2>r4 ^j^^ j^^g^ twofold 

integrals, which can be written as 



/ dT2 / dri e 
Jo Jo 

(■0 rli 

= dn d7-2e'^i("i-^3)+*^=("=-^*). (A14) 

Jo Jti 

By combining the last two integrals, (A13)and (A14), we 
obtain the law of conservation for the frequencies, 

^Tl>T3>T2>T4 _|_ ^T2>T4>ri>T3 
^Xi,X3 ^X2,X4 

The same procedure can be used for the other 7 groups 



,iTi(0i-03) +1x2(02 — 04) 



of integrals, which finally leads to 

^2(xi, if2i; X2, i02ix3, iOa; X4, iCti) 

= {^xi,X3 ^X2,X4 ^01,03 ^02, 

+ ^Xi,x4 ^X2,X3 fci,a4 ^2,03} 

(2^c)' 



This equation can be rewritten in the form. 



(A16) 



4)2 (xi , zfii ; X2 , iri2 |x3 , ifls; X4 , irJ4) 

= (l){-x.i , iOi|x3, i03) ^(x2, if22|x4, if24) 

+ (?i(xi,if2i|x4,i04) (^(X2,if22|x3,i03), (A17) 

0(xi,if2i|x3,iO3) = 4>{i9.i)l35ni,Q.3 '5xi,x3- (A18) 

From these equations we finally obtain Eq. (36). In sim- 
ilar manner we calculated the Fourier representation of 
the two-cloud function 932, which leads to Eq. (37). 
APPENDIX B: DETERMINATION OF THE 
GREEN'S FUNCTION r<,,<,,(iw) 

The evaluation of T^^^' {i'^) requires knowledge of the 
following two- and three-particle Green's functions: 

G°[a,T;a',r'|Ti] = (T a,(T) (r') n(Ti))o, (Bla) 

Gl[a,T;ai,T'\Ti,T2] = {T acr{T)aa'{T')n{Ti)n{T2))o, 

(Bib) 

where the statistical averaging is determined by the local 
part of the Hubbard Hamiltonian. In the calculation it is 
necessary to switch over to Hubbard transfer operators, 
X™", where the indices m and n denote ionic quantum 
states, m = 0, ±1, 2: 



Oct 



aO 



(B2) 



We omit calculational details and list only the resulting 

expressions. 



G0[c7,r;a',r'|n] = ^{^(r - r')^?(r' - n) e-'^^-+(---')^ 

+ 6{t - n)9in - t') [e-/5So+(r-r')A ^ ^^-fiE^{r-r')A^ 

- 9{t' - t)9{t - n) [e-'5^"+(^-^')^ + 2e-'3^^+(^-^')^] - 0(t' - n)e{T, - r) e-/3S2+(r-r')A 

+ e{n - t)6{t - t') e-/3^V+(r-r')A _ g^^^ _ ^y^^, _ ^g-/3£<, + (r-r')A ^ ^^-i3E, + ir-r')A^ ^ (B3) 

G°[a, r; a', r'ln, ra] = ^-^{[e{T - T')e{T' - n)e{n - T2) + 9{t - t')9{t' - T2)9{t2 - n)] e'/'^'^C-OA 
+ {9(t - n)9in - t')9{t' ~ T2) + 9{t - T2)9{t2 - t')9{t' - n)] e-/3^V+(--')A 

+ [9{t - Ti)9{Ti - T2)9{t2 - t') + 9{t - T2)9{t2 - Ti)9{n - t')] [e-/3^o+(r-r')A ^ ^^-0E^(r-T')A^ 
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— 




- t)6{t - 


Tl)9{Ti - 


T2) + e{T' - 


T)e{T- 


r2)e{T2 - 




— 


[6(t' 


- Ti)eiTi 


- T2)e{T2 


- r) + e{T' - 


- T2)e{T2 


- Ti)e{Ti 


- r)] e-''-^2+(T-T')A 


— 




- n)e{n 


-r)e{T- 


T2) + e{T' - 


T2)%2 


- T)e{T - 


n)] 2e-^^^+(---')^ 


+ 




-r)e{T- 


T')e{T' - 


T2)+6{t2- 


r)e{T- 


T')e{T' - 




+ 




-T)e{T- 


T2)0{t2 - 


t') + e{T2 - 


r)6»(r - 


Ti)e{Ti - 




+ 




- T2)0{t2 


- T)e{T - 


t) + e{T2 - 


ti)9{ti 


- T)e{T - 






[6»(ri 


- r')e{T' 


- T)e{T - 


T2) + 0{t2 — 


r'Wir' - 


- T)e{T - 


Tl)] [e~'^^'' + (^~^ >^ + 4e-P^2+(r-r )Aj 






-T')e{r' 


- T2)6{t2 


- r) + e{T2 - 


-r')e{T' 


- n)ein 








- T2)%2 


- T')e{T' 


- r) + e{T2 - 




- T')e{T' 


- r)] [e-/3E<^ + (^-^')A + 4g-/3B2+(r-r')Aj | ^ (34) 



where 

A = Eo-E^ = ^l, A = E^-E2=i^-U, (B5) 

and 0(a;) is here the step function. Integration over ti 
and T2 yields: 

Mr-r') = / rfTi(Ta<,(r)a<,.(T')n(ri))o 
= %^(0(t - r') (r - /) e-Z'^'+C— 
+ 9{t - r') (/? + t - r') e-/5^-+(^-')S 



-%'-r) (/3 + T-r') 
-6'(r-r') (2/3 + r-T')e 



_,/N„-/3B„ + (r-r')A 



-/3B2+(r-T 



')^}, (B6) 



,(r-T')= / rfn / dT2{Ta^{T)aa'{T')n{Ti)n{T2))o 
Jo Jo 

^f{e{r-r'){r-rr 



^2 g-/3Bo+(r-r')A 



/n2 /3£V+(r-T')A 



+ ^^(r - r') (/3 + T - t') 

- e{T' - r) (2/3 + r - r')' e-''^^+("-"')^}. (B7) 
The Fourier representation of these functions, 



Mr-r') = ^Y. 



IpniioJn), (B8) 



are given by 



tp2{ioj) 



■/3 



A(iw) A(ia;) 



(B9) 



Sa,a' [ 2(e 

^0 I 
2(e-^^ + e 



A3(iw) 



+ 



2/3 e- 



A^(ia;) 



2/3 e- 



4/3 e 



A2(itj) 

/32 e-('^' 



A (iw) 



A (iw) 



/32e- 



4/32 e- 



-/3B2 



A('iu;) X{iuj) A(zcj) 

This allows to evaluate the Tcr,a' {r — r') function, 



(BIO) 



r„,a'{r - t') = Mr - r') - 2{n)of3Mr - r') 

+ P'Gl,,{T-T'){{n\-2{n)l), (Bll) 

and to obtain its Fourier representation, 

<5<7,(7 



r<r,<T'(«'^) = 



2(e-/3^^o+e-^^-) 



Zo 
2(e-/5^ 



+ e 



A3(ia;) 

-/3E2) 



A^(ia;) 



+ 
+ 



2/3 



A2(ia;) 
2/3 



(n)o(e" 



+ e 



_2 [e-'^^ + 2e-^^= 
A (iw) 



-(n)o(e-^^ + e-'5^^)] 
A(zw) 

+ ((n2)o-2(n)2)(e-/5^«+e-'^^-)] 
+ =^[-(l-2(n)o)e-'^^ 
+ -(l-(n)o)4e-'5^^ 
+ ((n2)o-2(2n)^) (e'^^^ + e"^^^)] |. 



(B12) 
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For the case of half fiUing we have, 

Zo = 2(1 + e^**), X{ico) = ioj + n 
\{iu)) = iu) — fi, G°(ia;) = 



(B13) 



(ia;)2-^2 [(jw)2_^2]2 

/3m 



I 



V'3(«t^) 



(B14) 



2(1 + 6/3/^) 
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4/?[(ia;)2 + e^^ 4/3[(iw)2 + + 2iiOiJ.] 



[(ZC^)2 - m2]2 [(iw)2 - ^2]2 

2ia;/32e^'^^ 4/32(,^' + //) 
(ia;)2 — /i2 (ia;)2 — /i2 



(B15) 



The corresponding analytical continuation iui ^ E + 
iO'^ for E = gives to, 



MO) 



J3_^ 

(5,,,. 2/3/ ^ /3m 



V2(0) = ^(l + 3^), (B16a) 



M^ 



(B16b) 
(B16c) 



G°(0 = 0. 

In this case the F function is of particular simple form, 
r40) = ^(l + ^)[l-(n)o]. (B17) 

APPENDIX C: THE CRITICAL TEMPERATURE 

With help of the Poisson summation formula, 

i ^ fiiun) = ^ Jjiz tanh {^Pz)f{z), (CI) 

where C is the usual counterclockwise contour of the 

imaginary axis, wc obtain for case that the parameters 
in Eq. 123 obey 7^ > {c^ - Sfj^, 



1 



X Im 



/3 ^ Ql (iw) ^4^4 _ (^2 _ ^2)2 

/ tanh ^/3(ai + 102) 
V ai + ia2 

- rf2 tanh i/3c 1 
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2C74 2C74 



X Re 



tanh ^/3(ai + 10:2) 



c2-d2 tanhi/3(i 1 
^ ~ 2^^' ^ 2dy 



X Re 



tanh \d(Qi\ + 202) 



(C2a) 



(C2b) 



(C2c) 



with 



and 



(C3) 



Re 



Im 



tanh \(i{oL\ + 102) 
ai + ioi2 

ai sinh /3q!i + a2 sin /3q;2 
(a2 -|- q;2 -j [(.Qyjj ^jQ,^ _|_ cos /3a2 ' 

tanh ^/3(q!i + 20:2) 

ai + 7'Q'2 

«2 siiih .ini + oi riiu Ja2 



(C4) 



— • (C5) 

{af + Q!2)[cosh^Q!i + cos /3a2 ' 

Inserting the results in Eq. (C2) into Eq. (123) allows 
us to obtain expressions for the parameters 7711, 7722 and 
rii2. Furthermore, this allows to decompose the Tc equa- 
tion with the help of the following abbreviations, 

/ £/2 (i + e/3^)(e^^ + e^(2,.-c/)-) X 

"V i'2^^-ur zl ) 

X (?7ll?722 - ?7i2) 

_/ \ (i + e^'')(e^'' + e^(2^-^)) 

X (t,r^.c)2 (c'-f')' f(' tanhl/3c ^ tanh 1 /3d 



478 [ \ c 
/tanh 5/3(ai + 10:2) \ tanh i /3c tanh ^ /3c 

X xxC I I 

V ai+ ia2 J c c 

_ / tanh|/3(ai+iQ2) Y 27^ - (c^ - rf2)2 
V ai + ias y ^474 - (c2 - rf2)2 

/ tanh iB((y-i + i(y.-^)\ 

X Im 



/ tanh ^/3(ai + ia2)\ f tanh i/3rf tanh i/3c 
\ 0:1+ ia2 J \ d c 



(C6) 



and 

B={1 + I — 7711 



2^1- U J Zo 

U \ e^^ + e'3(2^-c/) 
+ I 1 - 77 I 7i ??22 



2ii-U J Zo 
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-l-Im 



Zo{2n-U) 



X IU + 



/ tanh ^/3(q:i + ia2) 
\ ai + ia2 

(2M-C/)(l-e'5(2M-^)) 



tanh i/3c (1 + e^'') 



(C7) 
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where the value of Ii is determined by Eq. (C2a) . Equa- 
tion (105), which determines T^, can then be written as 

1 + Bf, + Af^ = 0, (C8) 

where the quantity is given by Eq. (96). Equation 
(C7) would allow to investigate in detail the interplay 
of the different parameters and renormalized quantities 
obtained after eliminating the electron-phonon interac- 
tion by the Lang-Firsov transformation. The influence 
of these parameters on Tc can however only be obtained 
by numerical work, which has yet to be undertaken. The 
simplified discussion in the main text shows that Tc is 
proportional to the typical energy scale involved, which 
is the renormalized chemical potential and not the bare 
quantities of the original model. This means that due 
to the Lang-Firsov transformation the proportionality of 
Tc to a typical phonon frequency as in BCS theory or in 
the Eliashberg formulation is lost. This is an interesting 
observation but must be tested numerically by analyz- 



ing the more complex expressions of this paper. In order 
to find out whether the proportionality Tc cx i.e., 
Tc proportional to a renomalized electronic energy scale 
and not proportional to an average phononic frequency, 
(wk), is an intrinsic property of the model, we will per- 
form investigations in infinite dimensions. This is nowa- 
days much at debate by using other approaches like the 
dynamical mean field theory, which allows to evaluate 
more or less accurate materials properties if combined 
with density functional theory'^^"'^^. We believe that our 
approach allows for similar accurate description of ma- 
terials including superconducting properties. This is left 
for future work. 

Finally we would like to stress that the present ap- 
proach, to start from the exact local (atomic) description 
and to take into account the properties of the correlated 
tight-binding electrons on a lattice by perturbation the- 
ory in the transfer integral, has conceptually some ad- 
vantages compared to other theories. 
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